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Motivated by the ATLAS and CMS discovery of a Higgs-like boson with a mass around 125 GeV, 
and by the need of explaining neutrino masses, we analyse the three canonical SUSY versions of 
the seesaw mechanism (type I, II and III) with CMSSM boundary conditions. In type II and III 
cases, SUSY particles are lighter than in the CMSSM (or the constrained type I seesaw), for the 
same set of input parameters at the universality scale. Thus, to explain m h o ~ 125 GeV at low 
energies, one is forced into regions of parameter space with very large values of mo, M1/2 or A)- We 
compare the squark and gluino masses allowed by the ATLAS and CMS ranges for m h o (extracted 
from the 2011-2012 data), and discuss the possibility of distinguishing seesaw models in view of 
future results on SUSY searches. In particular, we briefly comment on the discovery potential of 
LHC upgrades, for squark/gluino mass ranges required by present Higgs mass constraints. We also 
take into account the MEG upper limit on Br(/z — > e~f) and show that, in some cases, this may help 
to restrict the SUSY parameter space, as well as to set complementary limits on the seesaw scale. 
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I. INTRODUCTION 



With the data accumulated in 2011 and 2012, both the CERN ATLAS and CMS col- 
laborations have recently claimed the discovery of a new particle that resembles very much 
the long-awaited Higgs boson. The mass of this new state, measured in good accordance in 
different decay channels, is in the ballpark of m^o ~ (123 — 127) GeV. While the overall 
significance in the 2011 data was only 2.2cr in ATLAS |lj and 2.1a in CMS j^], with the 
2012 update both experiments increased their statistical significances to the 5cr discovery 
threshold jsl, 0|. Especially noteworthy is that both ATLAS and CMS observe an excess of 
events in the 77 and ZZ decay channels with an invariant mass which differs by roughly 2 
GeV, i.e. consistent at the la level. Complementary evidence has been reported by the CDF 
and DO experiments at the Tevatron. These collaborations have recently released updated 
combined results on searches for the Higgs boson 5j, finding a ~ 3a statistical significance 
in the bb decay channel. 

Given that supersymmetry (SUSY) has been the most popular paradigm for physics 
beyond the standard model (SM) in the last decades, the recent LHC results have triggered 
the expected flurry of theoretical activity dedicated to the study of how a relatively heavy 
Higgs constrains the supersymmetric parameter space . The general consensus is that 

a lightest Higgs boson with a mass of m h o ~ 125 GeV is uncomfortably heavy for minimal 
SUSY. Here, by minimal SUSY we mean a supersymmetric model with no new superfields 
and no new interactions, gauged or non-renormalizable, at the electroweak scale. In this 
ramework, the hefty Higgs requires either multi-TeV scalar tops or very large stop mixing 



16, 31 



[32]. In the latter case, the lightest stop could still be relatively light, say 
m t 1 ~ 500 GeV (43J. For a constrained minimal supersymmetric standard model (CMSSM) 
with universal boundary conditions at a high scale, such a spectrum requires that at least one 
of the three basic parameters M1/2, m or A takes a minimum value of several TeV (7, 10, 31 



In addition, it has been found that a moderately large tan f3 may be helpful to increase the 
Higgs boson mass 0, EHj • 

The naturalness problem of the MSSM with a ~ 125 GeV Higgs mass has revived the 
discussion around non-minimal supersymmetric extensions of the standard model. In par- 
ticular, the recent LHC data has been scrutinized in the context of SUSY models with new 
like the NMSSM) 15J, |22j, |26j, |35|, |39|, extended gauge models with additional new 



F-terms 
D-terms 



44H49|, heavy-SUSY scenarios like Split SUSY 



50- 5j "natural SUSY" 



33, 53-55J 



and high-scale SUSY 56|, |57[) or "effective" SUSY, i.e. SUSY with new non-renormalisable 



operators 58M61]. among others. In this work, we will follow an alternative approach and 
assume SUSY is realized minimally. We explore the consequences of the LHC Higgs search 
data on the CMSSM parameter space and the SUSY spectrum, from a viewpoint similar 
to that taken in MSSM-dedicated studies like, for instance, the one of Ref. hOll. However, 
our analysis differs from these by considering that a seesaw mechanism for neutrino mass 
generation is implemented in the MSSM. Our motivation lies in the fact R-parity conserving 
MSSM (with or without CMSSM boundary conditions) does not provide an explanation for 
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the observed neutrinos masses and, thus, is not complete. 

From the theoretical point of view, implementing the seesaw mechanism in the (super- 
symmetric) SM seems to be the simplest (and most motivated) solution to the neutrino 
mass problem. With renormalizable interactions only, there are three tree-level realizations 

These 



of the seesaw mechanism |62| usually called type I |63j-|65j, II |66H72f and III 
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three variations differ from each other by the nature of their seesaw messengers. Namely, 
in type I an effective neutrino mass operator arises from the decoupling of heavy neutrino 
singlets, while in type II one integrates out a heavy SU(2) scalar triplet with hypercharge 
two. Instead, in the type III seesaw neutrino masses are generated through the tree-level 
exchange of SU(2) fermionic triplets of zero hypercharge. If in type II and III one extends 
the MSSM by just adding the superfields required to generate neutrino masses, then one 
of the most appealing properties of the MSSM is lost: gauge coupling unification. This 
stems from the fact that both the scalar and fermionic triplets belong to incomplete SU(5) 
representations. Unification can be easily restored by embedding those states in full SU(5) 
multiplets like 15-plets in the case of type II 74| or 24-plets 75] in the case of type III. 
Note that, in addition to the SU(2) triplet, the 24 of SU(5) contains a singlet which also 
contributes to the effective neutrino mass operator and, thus, the decoupling of the 24-plet 
leads to an admixture of type I and type III seesaws. 

The main purpose of this work is to investigate whether imposing a Higgs mass around 
125 GeV allows to some extent to differentiate the CMSSM from the constrained SUSY 
seesaws and also whether type II and III seesaws are distinguishable among themselves. We 
will complement this analysis by imposing the MEG constraint on the branching ratio of 
the radiative lepton flavour violating decay Br(/i — > e^y) < 2.4 x 10~ 12 
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The rest of this paper is organized as follows. We start by recalling the general features of 
the aforementioned SUSY seesaw models in Section [XT] and present some discussion related 
with lepton flavour violation (LFV) in Section IHIl Afterwards, we describe our numerical 
analysis and present its results in Sections llV Al and IIV Bl respectively. Our conclusions are 
drawn in Section IVl 



II. MODELS 



In the following we will briefly describe the three types of SUSY seesaw mechanisms 
considered in this work and possible embedding in a grand-unified (GUT) model based on 
the SU(5) gauge group. We use standard notation for the MSSM superfields, namely L, Q 
and H u (Ha) denote the lepton, quark and hypercharge one (minus one) Higgs superfields, 
while the lepton and quark singlets are E c , D c and U c . The vacuum expectation values of 
H ud are denoted by v U)d /\/2 with tan/3 = v u jv d and v = v\ + v\ = 246 GeV. 
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A. Supersymmetric seesaw type I 



In the case of the supersymmetric type I seesaw, very heavy singlet superfields N c are 
added to the MSSM, yielding the following superpotential below the grand-unification scale 
M GUT : 



Wi 

Wmssm 



w. 



MSSM 



+ w v 



Y U U C QH U - Y d D c QH d - Y e E c LH d + ^H u H d 
Y V N C LH U + l -Wl R N c N c , 



(1) 
(2) 

(3) 



where SU(2)-invariant products are implicit. This model can be realized in an SU(5) frame- 
work taking the following SU(5) matter representations: 1 = N c , 5m = {D c , L} and 
10m = {Qi U c , E c }. At the effective level, a dimension five neutrino mass operator of the 
type LLH U H U originates from the decoupling of the heavy singlets, leading to an effective 
neutrino mass matrix given by the well-known seesaw formula 

m„ = -^Y^M^Y,, (4) 



after electroweak symmetry breaking (EWSB). Being complex symmetric, is diagonalized 
by a 3 x 3 unitary matrix U 

(5) 
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U T m u U . 

The lepton mixing matrix U can be parameterized in the standard form 



/ 



U 



C12C13 



S12C13 



0125235136'" 
i5 

\ s 12 s 23 ~ Ci2C23Si3e 



si 3 e 

C12C23 - s 12 s 2 3Si3e lS S23C13 
id 



is \ 



/ gtai/2 







V 








e i* 2 /2 Q 

1 



, (6) 



-S12C23 

AX 

"C12S23 — Si2C23Si3e C23C13 J 

with Cij = cos 8 \j and s^- = sin^jj. The angles 6 ) i 2 , ^13 and 6*23 are the solar, the reactor 
(or CHOOZ) and the atmospheric neutrino mixing angle, respectively, while 5 is the Dirac 
phase and 0:1,2 are Majorana phases. 

It is well known that the Dirac neutrino Yukawa couplings Y u can be defined in terms of 
the physical neutrino parameters, up to an orthogonal complex matrix R 
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Y v = ^-VMfiRVrivU 



(7) 



where m„ and Mr are diagonal matrices containing the light and heavy neutrino masses, 
respectively. It is worth noting that, in the special case of R = 1, the non-trivial flavour 
structure of Y u stems from the lepton mixing matrix U. 

B. Supersymmetric seesaw type II 

In the type II seesaw, neutrino mass generation is triggered by the tree-level exchange of 
scalar triplets. Its simplest SUSY version requires the addition of a vector- like pair of SU(2) 
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triplet superfields T and T of hypercharge Y = ±2. A natural way to implement the type 
II seesaw in a GUT scenario is to embed the triplets in a 15 and 15-plet of SU(5) which 



decompose under SU(3) x SU(2) x U(l) in the following way |74 | 

15 = S + T + Z , (8) 
S ~ (6, 1, -2/3), T ~ (1, 3, 1), Z ~ (3, 2, 1/6), 
with an obvious decomposition for the 15. The SU(5) invariant superpotential reads 

W = -=Y 15 5 15 5 + -=X 1 5 H 15 5 H + -^A 2 5 H l5 5^ + Y 5 10 5 5 H 
v2 \/2 v2 

+ Y 10 1010 5h + M 15 15T5 + M5 5h5^, (9) 

with 5h = {H C ,H U ) and 5# = (H c ,Hd). We do not go through the details of the SU(5) 
breaking as we take the above SU(5) realization only as a guideline to fix some of the 
boundary conditions at Mqut- Below Mqut, in the SU(5)-broken phase, the superpotential 
reads 

W n = Wmssm + -^(Y T LTL + Y S D C SD C ) + Y Z D C ZL 
v2 

+ -^(X.HjTHj + X 2 H U TH U ) + M T TT + M Z ZZ + M S SS . (10) 
v 2 

The dimension five effective neutrino mass originates now from the decoupling of the triplet 
states, leading to an effective neutrino mass matrix 

^2 



vl A; 



m u = — — — , (11) 
2 M T 

once electroweak symmetry is spontaneously broken. It is apparent from the above equation 
that the flavour structure of m u at low energies is the same as the one of the couplings Yt 
at the decoupling scale Mt (up to renormalization group effects which can be relevant in 



some special cases 78[). Consequently, Yy is diagonalized by the same matrix as m^, i.e. 

Y T = U T Y T U. (12) 

In short, if all neutrino eigenvalues, angles and phases were known, Yt would be fixed up 
to an overall constant which can be easily estimated to be 

M T in /0.05eV 



- 10 15 GeV . (13) 

A 2 V m„ J 

In principle, the remaining flavoured Yukawa couplings Y $ and Y z are not determined by 

any low-energy neutrino data. Still, they both induce LFV slepton mass terms, just as Yt 

does. Having the above SU(5) GUT model in mind, we impose the unification condition 

Yt = Y,s = Y z at Mqut in our numerical analysis presented in Section IIVI As for the 

heavy-state masses, the mass equality condition Mt = Mz = Mz = M15 imposed at the 

GUT scale is spoiled by the renormalization group (RG) running of the masses. Nevertheless, 

these effects are small and, therefore, gauge coupling unification is maintained. In view of 

this, for practical purposes we decouple the T, Z and S states at the common scale Mt(Mt), 

neglecting in this way threshold effects resulting from the small RG-induced splittings among 

the heavy masses. 
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C. Supersymmetric seesaw type III 



In the case of a type III seesaw model, neutrino masses are generated by the tree-level 
exchange of zero hypercharge fermions, usually denoted as E, belonging to the adjoint repre- 
sentation of SU(2). These states can be accommodated, for instance, in a 24-plet of SU(5). 
Above the SU(5) breaking scale, the relevant superpotential for our discussion is 

W = y/2Y 5 l M lQ M l H - ^Y 10 10 M 10 M 5^ + Y 24 5 H 24 A/ 5 A/ + 1m 24 24 m 24 m . (14) 

As in the type II case, we do not specify the Higgs sector responsible for the SU(5) breaking. 
The superpotential terms directly involved in neutrino mass generation are those containing 
the representations 24^, which decompose under SU(3) x SU(2) x U(l) as 

24 M = (1,1,0) + (8, 1, 0) + (1, 3, 0) + (3, 2, -5/6) + (3*, 2, 5/6) , (15) 

= n c + g + z + x + x. 

The fermionic components of (1, 1,0) and (1,3,0) have the same quantum numbers as N c 
(the type I heavy neutrino singlets) and E. Thus, one expects that, in general, the decoupling 
of the 24m components leads to an effective neutrino mass operator which contains both a 
type I and a type III seesaw contribution. In the SU(5) broken phase the superpotential is 

Wui = Wmssm + H u | Y S E - ^Y U N C | L + Y X H U XD C 

+-m r n c n c + -m g gg + -m^j: + m x xx . (16) 

2 2 2 

Once more, we impose the GUT scale boundary condition Y^ = Y v = Yx and = 
Mg = = Mjf. Integrating out the heavy fields, and after EWSB, the following effective 
neutrino mass matrix is generated: 



m,. 



f (^M^Y, + iY|M^Y s ). (17) 



As mentioned above, there are two contributions to neutrino masses stemming from the 
gauge singlets N c as well as from the SU(2) triplets E. In this case the extraction of 
the Yukawa couplings from low-energy parameters for a given high scale spectrum is more 
complicated than in the other two types of seesaw models. However, as we start from 
universal couplings and masses at Mqut-, we find that at the seesaw scale one still has 
~ Ms and Y v ~ Y^;. Consequently, one has 

* -vl ^ YljM^Yjj , (18) 

to a good approximation. This result allows us to use the same decomposition for Ye as 
the one discussed in section |H"A~| up to the overall factor 4/5 [see Eq. ([7])]. 
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III. LEPTON FLAVOUR VIOLATION IN THE (S)LEPTON SECTOR 



The search for LFV processes beyond neutrino oscillations has attracted a great deal of 
attention both from the experimental and theoretical communities. Rare decays like fi — y 
have been searched for decades, without any positive result. The most stringent constraint 



on this process comes from the MEG experiment [79| which, by analysing the data collected 



in 2009 and 2010 
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has set the new bound Br(/i — y ej) < 2.4 ■ 10 



-12 



The branching ratio (BR) for ^ — > Ipj can be generically written as 



80 



4S7r 3 (T 

Br(/, Xfi) = {\Al\ 2 + |A^| 2 ) Bx{k Iju^) . (19) 

The amplitudes and Ar depend on the specific physics framework and, in general, are 
generated at the 1-loop level. In our SUSY scenario, the dependence of those amplitudes on 
the LFV slepton soft masses is approximately given by 

Ag~^ , A«~^, (20) 

m SUSY m SUSY 

where m| and m| c are the doublet and singlet slepton soft mass matrices, respectively, and 
msuSY is a typical supersymmetric mass. In the derivation of these estimates one typically 
assumes that (a) chargino/neutralino masses are similar to slepton masses and (b) left-right 
flavour mixing induced by A-terms is negligible^. 

Assuming universal boundary conditions for the soft SUSY-breaking terms at the GUT 
scale, and considering only the leading-log approximation for the LFV slepton masses and 
trilinear terms induced through RG running, one obtains: 

/ 2\ a k 



{3ml + <) ( Y 1 LY ^. > ( 21 ) 
( A *h * A) (Y e YiLY fc ) _ _ , (22) 



Llv ~ 8tt 2 



16tt 2 

for % 7^ j. In the basis where Y e is diagonal, L mn = \ii(Mgut / M n )S mn and is the Yukawa 
coupling of the type-k seesaw (k = I, II, III) with Y^ = (Y^Y^Y^), given at Mqut- 
Taking into account the renormalisation group equations (RGEs) for m| and A e we obtain 

ai = 1 , an = 6 and am = 9/5 . (23) 

Note, that in case of the type II seesaw the matrix L is proportional to the identity and 
thus can be factored out. All models considered here have in common that they predict 
negligible flavour violation for the RH sleptons 

(mlOy * 0. (24) 



1 This assumption is not valid when large values of \Ao\ are considered. Nevertheless, the above estimates 
can be still used to illustrate the dependence of the BRs on the low-energy neutrino parameters. 
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Although not very accurate, the above approximations allow to estimate the LFV slepton 
masses and A-terms within different seesaw frameworks. The BRs for rare lepton decays 
li —7- Ijj are roughly given by 

Km 2 -)- -I 2 

Br (li -»• L- 7 ) oc a 3 ml R L 3 tan 2 f3. (25) 

rn SUSY 

For distinct seesaw scenarios, and a given set of high-scale parameters, the above BRs change 
due to the different (m^) 2 - and the distorted mass spectrum (which differs from the pure 
CMSSM one). The most important parameter turns out to be the seesaw scale due to its 
influence on the size of the Yukawas. The higher the seesaw scale is, the larger are the 
Yukawa couplings and, consequently, the LFV rates. In case of the type II seesaw, the 
coupling A2 plays a crucial role, as seen in Eq. (fT3|) . Small values of this parameter lead to 
large Yy Yukawa couplings and high LFV rates. 

Finally, we would like to comment on the influence of the R matrix on LFV decay rates. 
As shown in Eq. (J7|), the ~Y U Yukawa couplings for type I seesaw are proportional to R and, 
thus, different choices of this matrix lead to different off-diagonal entries in the soft squared 
mass terms (which in turn changes the LFV rates). Similarly, the type III Yukawa couplings, 
Ye, follow an analogous equation and, consequently, also change with R. This additional 
freedom can be used to cancel some (y^ISY/A combinations, in particular the one with 

(i, j) = (/i, e) (77j|. This allows for large LFV effects in the r— e and r—fi sectors while having 
negligible \x — e transitions. In the following, we will disregard this possibilitjj 2 ]. Therefore, 
implications on M$s drawn from \x — > e'-f considerations can be regarded as approximate 
lower bound^E 



IV. NUMERICAL ANALYSIS AND RESULTS 



A. Setup 



S3 



Our numerical results have been obtained with SPheno 
parameters, as well as the usual universal soft terms at the GUT scale 

m , Mi/2, A), tan/3, sign(/x), 



84j. Taking as input the SM 

(26) 
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SPheno computes the resulting SUSY spectrum by means of complete 2-loop RGEs 
properly adapted for every model. This includes the pure CMSSM and the three seesaw 
variants studied in this work. At the SUSY scale, the /1 parameter is obtained including the 

2 In fact, we will always consider real parameters and degenerate spectra for the right-handed (RH) neutrinos 
(in type I) and for the SU(2) fermion triplets (in type III). In such scenarios the R matrix is physically 
irrelevant, since it drops out in the computation of (St^LY/.) 77|. For a discussion on the effects of 

0,|82|. 

1. Larger Mss are in 

principle possible if one tunes the R matrix in order to have a cancellation in the fx — e sector. However, 
one cannot find R matrices that allow to go to much smaller Mss scales. 



considering complex parameters we address the reader to, e.g. Refs 
Given Br(^t — > 67) and mo ~ msusY one can determine Mss assuming R 
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most relevant 2-loop corrections [88j and complete 1-loop corrections to all sparticle masses 



are implemented 89j. These calculations follow the DR renormalization scheme. 

In case of the Higgs boson mass, the aforementioned 1-loop corrections are supplemented 
by the most relevant 0[a s (a t + a&) + (a t + a^) 2 + a T ab + a 2 .] 2-loop contributions 



88, 90-9 



For a detailed study of the SPheno results for the Higgs boson mass and a comparison to other 



popular numerical codes we refer to 94]]. We have checked that our results agree, within 
the usual 2 — 3 GeV theoretical uncertainty, with the results given by FeynHiggs 95). This 
code uses an on-shell renormalization scheme and therefore small differences are expected 
on theoretical grounds. In particular, larger differences are found for very large Higgs boson 
masses, m^o ~ 129—130 GeV, a region where numerical computations are no longer accurate. 
We also found good agreement between our results and those presented in some recent works 
devoted to the study of the Higgs mass in the MSSM 0M |3l|, |41|. Although the theoretical 
error is always present, and the exact numbers might differ in some cases, the general 
behaviour and the dependence on the SUSY parameters are correctly reproduced. We have 
decided not to compute the SUSY spectrum at a fixed scale Q = 1 TeV, as suggested by 
the SPA conventions (96j), since that is known to give a poor accuracy in the determination 
of the Higgs boson mass for scenarios with very large values of Aq or with multi-TeV stop 
masses. Instead, we compute the SUSY spectrum at the geometric average of the two stop 
M s 



masses 2 , i.e. 



Although we evaluate the Higgs mass numerically taking into account the higher-order 
corrections enumerated above, we find it useful to recall that the leading 1-loop corrections 
to the Higgs mass for m oderate values of tan (3 and large Higgs pseudoscalar mass tua, are 
approximately given by 



97- 



I5J 



2 2 n 3?77y 

m h o ~ m z cos p + - 

4n z v z 



In 



2 

mf 



+ 



Ml 



xl 

12M 2 



X t = A t - ncotp, (27) 



where \x is the Higssino mass parameters, A t is the top trilinear term at low-energy and X t is 
the mixing parameter in the stop sector. Obviously, the above approximation is not always 
very accurate. In any case, we will only use it to understand the behaviour of the Higgs mass 
with some of the input parameters of the seesaw models discussed in the previous section. 

In the following we will present and discuss our numerical results. Notice that we will 
loosely talk about "the seesaw scale", Mss, when referring to the mass of the seesaw medi- 
ators, i.e. the right-handed neutrino mass, Mr, in case of seesaw type I, the Y = 2 triplet 
mass, M15 (or Mr), for type II or the mass of the Y = triplet, M24 (or Ms), for type III. 
Our assumptions regarding the input parameters for each of the seesaw models are: 

Type I: We consider the general case of 3 degenerate RH neutrinos with mass Mr. In 
the flavour sector, we fix R = 1 [see Eq. (0)]. This choice does not have a significant impact 
on the Higgs mass since, as already pointed out, the effect of the Yukawa couplings on m h o 
is marginal. 

Type II: Apart from the unification conditions for the Yukawa couplings and masses of 
the different 15-plet components mentioned in Section [II Bl we will use in most of the cases 
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the values Xi^Mgut) = 0.5 for the superpotential couplings of the triplets with the Higgs 
superfields. Later, we will comment on how relaxing this condition affects the Higgs and 
squark masses. 

Type III: We will always assume the existence of three copies of 24-plets, with an 
approximately degenerate mass Ms. Alternatively, one could also explain neutrino data 
with two degenerate 24's or with three, being one "light" and the other two close to the 
GUT scale. The first of these options leads to results somewhere between those shown for 
type II and type III with three degenerate 24's, while the latter has m^o somewhere between 
type II and type I. Since nothing qualitatively new results from these cases, we will not 
discuss them in detail. As in the type I case we assume R = 1. 

For all our numerical cases the values of the low-energy neutrino parameters (mixing 
angles and mass-squared differences) coi ncide wit h the best-fit values provided by global 
analysis of all neutrino oscillation data 101 103| . To simplify our analysis, we consider 



all couplings and mass parameters to be real and for tan (3 we take the reference value 
tan/3 = 20. For other values of tan/3, our CMSSM results agree quite well with those 



discussed, for example, in [10|. We have scanned the parameters mo and Mu2 in the range 
of [0, 10] TeV. As for A , we have taken values in the interval [—5, 5] TeV, although we will 
mainly concentrate on the two extreme cases with A = TeV and A = — 5 TeV. For other 
choices of Aq, the results always lie between the extreme ones. 



Current bou nds on squark and gluino masses in CMSSM-like setups from ATLAS [10 



and CMS 105| already exclude m g = nig ~ 1.4 TeV and m g ~ (800 — 900) GeV for very 
heavy squarks. In our plots we therefore show only that part of parameter space where m g 
and rriq are larger than 1 TeV. 



B. Results 



It is well known that adding seesaw mediators with masses between the SUSY and GUT 
scales changes the RG running of gauge couplings. As a result, the RG flow of a ll Yukaw a 
couplings and mass parameters is modified with respect to the CMSSM case 106M110 . 
In the case of type II and III seesaws, the increase in the valu e of the common gauge 



106 



107( | . Therefore, one expects 



coupling o>(Mgut) leads, in general, to lighter sparticles 
the Higgs mass to be sensitive to the parameters characterising each seesaw model, namely 
the mass M$s and possible couplings with the Higgs and/or lepton sectors of the MSSM. 
Consequently, the reconstruction of the SUSY-breaking parameters at the universality scale 
Mqut from low-energy mass measurements will be very sensitive to the presence of new 
fields at intermediate scales. 

As an example, in Fig. [T] we show the behaviour of the stop masses m~ tl 2 and mixing 
parameter X t , as well as the mass of the lightest Higgs, h°, as a function of the seesaw 
scale Mssi for seesaw type I (blue), type II (red) and type III (black) taking m = 0.5 TeV 



and M- 



1/2 



1.5 TeV. The results are shown for two values of the common trilinear term 



10 



3.0 




M (GeV) M ss ( GeV ) 



FIG. 1. Variation of the scalar top masses (top panels), the ratio \Xt/Mg\ and the mass of the 
lightest Higgs h° (bottom panels) with M$s, for a particular point in the SUSY parameter space 
with mo = 0.5 TeV, Min = 1.5 TeV. The blue, red and black lines correspond to seesaw type I, 
II and III, respectively. The results are shown for Aq = TeV (dashed lines) and A$ = —3 TeV 
(solid lines). For values of the seesaw scale larger than roughly 10 15 GeV no solutions consistent 
with observed neutrino data can be found. Also, for M$s % (few) 10 9 [10 13 ] GeV gauge couplings 
become non-perturbative below the GUT scale in case of seesaw type II [type III]. See also text. 

at the GUT scale, namely Aq = —3 TeV (solid lines) and Aq = TeV (dashed lines). For 
values of the seesaw scale larger than roughly 10 15 GeV no solutions consistent with observed 
neutrino data can be found, while for values of the seesaw scale below approximately (few) 
10 9 (10 13 ) GeV gauge couplings become non-perturbative below the GUT scale in case of 
seesaw type II (type III). This explains why no results are shown for lower values of M$s in 
those cases. 

The first immediate (and expected) conclusion that one can infer from the results pre- 
sented in this figure is that there is essentially no dependence of m h o on M R in case of type 
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I (bottom-right panel in Fig. [I]). This is due to the fact that sparticle masses do not change 
with Mr, as can be seen for the particular cases of mf (top panels) and neither stop mixing 
does. Due to the singlet nature of the RH neutrinos in the type I seesaw, and to the fact 
that they only have Yukawa couplings Y u with the lepton and Higgs doublet superfields [see 
Eq. (DO)] the soft SUSY breaking MSSM parameters affected at the 1-loop level are m 2 Hu , m|, 
A e and A u . Still, even those show only very mild departures from their CMSSM values. All 
other soft masses change only at the 2-loop level. 

The neutrino Yukawas Y " v required to fit neutrino data depend on Mr and are 0(1) for 
Mr ~ 10 15 GeV. The results of all plots in Fig. [T]show that, even for such large Yukawas, 
the changes of the SUSY spectrum are relatively small, due to the short RG running from 
the GUT scale to Mr. For smaller values of Mr no traces of the seesaw remain in the SUSY 



spectrum The only important consequence of changing Mr is the strong effect of this 

scale on the LFV entries of m| which control the rates of LFV processes like \x — > ey (see 
Section [TTTI and the discussion below). Moreover, as the results of Fig. [1] show, changing 
the value of Aq from zero to -3 TeV, shifts down the stop masses (middle left panel) due to 
the term proportional to A\ in the RGE of (m| c )33. On the other hand, the magnitude of 
the stop mixing parameter X t increases as a consequence of the fact that, at low energies, 
\A t \ is larger for Aq = —3 TeV than for Aq = TeV. Therefore, \X t /Ms\ increases when 
going from vanishing Aq to Aq = —3 TeV, resulting in an increase of the Higgs mass by 
approximately 3 GeV. Of course, this feature is also present in the CMSSM and is by no 
means related with the presence of the heavy neutrino singlets. 

The situation changes when one turns to the type II and type III seesaws. In these 
cases, for a fixed choice of CMSSM parameters, one usually finds that stop masses become 
smaller when lowering the seesaw scale, M$s- At the same time, the stop mixing angle can 
increase. However, this increase is practically never sufficient to compensate for the smaller 
stop masses. Thus, in general, m h o decreases with decreasing M55 for both type II and 
type III. As Fig. [1] shows this decrease depends also on A , with changes in m h o being much 
smoother (and smaller) for Aq = TeV than for A = —3 TeV. For the lowest values of 
Mss possible, m^o can be even lighter for Aq = —3 TeV than for Aq = TeV. This is due 
to the rather strong dependence of the stop masses on Aq. All these features, discussed here 
for a special CMSSM point, are qualitatively valid for rather larger ranges on the CMSSM 
parameter space, as we will discuss next. 

In Fig. [2] we show examples of squark mass, Higgs mass and Br(/x — > e'y) contours in 
the plane (r^M^) for CMSSM plus seesaw type I, taking two extreme values of the 
seesaw scale Mr, namely 10 9 GeV (top) and 10 14 GeV (bottom); as well as two values of A : 
Aq = TeV (left) and Aq = —5 TeV (right). The hatched regions on the left lead to a charged 
LSP and, thus, are not acceptable due to cosmological constraints (charged dark matter). In 
the grey regions, either EWSB is not possible in a consistent way or gluino/squark masses 
are below 1 TeV. The solid lines show contours of m h o at 124, 125.3 (central CMS value), 
126.5 (central ATLAS value) and 127 GeV, to reflect the currently favoured region of m^. 
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FIG. 2. Examples of squark mass, Higgs mass and Bi(p — > cy) contours in the plane (mo, M1/2) 
for CMSSM plus seesaw type I for two values of the seesaw scale M R : M R = 10 9 GeV (top) and 
Mr = 10 14 GeV (bottom); as well as two values of A : A = OTeV (left) and A = -5 TeV 
(right). Br(/i — > e^) is orders of magnitude below the expected experimental sensitivity in case of 
Mr = 10 9 GeV and, therefore, contours are not shown (for a discussion see text). 



The green dashed lines correspond to constant average squark masses, denned as 

mi + mi 

m _ = dL _ dR ^ (2g) 

while the (black) dash-dotted lines refer to the contours Br(/i — > c-f) = 2.4 x 10 -12 (MEG) 
and 10~ 13 . The yellow region corresponds to values of m^o in the CMS interval 125.3 ± 
0.6 GeV. 

For m h o = 125 GeV we find squark masses in the range of typically 5 TeV for A Q = TeV 
and as low as 2 TeV for A = —5 TeV. However, considerably larger squark masses, 
(9(10) TeV, can be found in the CMS preferred window of Higgs mass. This is consistent with 
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the findings of previous works on the CMSSM @,E1 

l3l| and in agreement with expectations. 
If this scenario is indeed realized in nature one expects to observe squarks at the LHC with 
y/s = 14 TeV (LHC14) only for the largest values of Aq. Still, even for Aq = —5 TeV large 
parts of the allowed parameter space in squark and gluino masses will remain unexplored 
by LHC14. 

Fig.[2]shows also contours of Br(/x — > e^y) assuming (a) a degenerate RH neutrino spectrum 
with R = 1 (see Section III Ap and (b) low-energy neut rino s fitted with a normal hierarchy 
spectrum and mixing angles within the allowed range 10l| . It is well-known tha t differen t 
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choices of #13 can lead to values of Br(/x — > ej) dif ferin g by a considerable factor 
However, we fix 13 according to the results of 10 1| . where a global fit to all available 
experimental data gives a best-fit value of sin 2 6*13 = 0.026 in case of a neutrino spectrum 
with normal hierarchy. We note that with degenerate RH neutrinos and R = 1, a complete 
canc ellation of Br(/i — > ej) is no longer possible within the 3cr allowed range of sin 2 # 13 



101] . In case of Mr = 10 9 GeV, Br(/x —> ej) is orders of magnitude below the expected 
experimental sensitivity and, thus, contours are not shown. However, if Mr = 10 14 GeV, 
Br (/1 — > ej) is well within the current expected sensitivity of MEG. Therefore, for the 
CMSSM with a seesaw type I and m^o = 125, MEG already provides an upper limit on Mr 
of the order of 10 14 GeV, despite the fact that sleptons in the CMSSM are relatively heavy 
in the allowed parameter space. We notice that the constraints from Br(/i — > ej) are, in 
general, more stringent for large values of Aq. 

In Fig. [3] we show the results in the plane (mo, M1/2) for type II seesaw with Mt = 
10 9 , 10 14 GeV and Aq = 0, -5 TeV. When M T = 10 14 GeV, the results for m h o and the 
squark masses are very similar to the CMSSM ones, although some small shifts are visible 
upon closer inspection (see also below). On the other hand, the cont ours f or B r(/i — > ej) are 
different from those in Fig.[2j This is in agreement with expectations 106l . ll 14] . since in type 
I neutrino masses scale as the square of the Yukawa couplings whereas in type II neutrino 
masses are linearly proportional to Yt [see Eqs. (jl]) and (flTT)]. while the RG running of the 
LFV soft masses depends quadratically on the Yukawas in both cases. Note, that in Fig. [3] 
we have used A2 = 0.5. A value of A2 = 1 would lead to smaller values of Br(/i — > ej) 
by (roughly) a factor of four. Much larger values of A2 are not allowed, if the theory is to 
remain perturbative up to the GUT scale. 

For Mt = 10 9 GeV, on the other hand, the results look drastically different. All Higgs 
(and squark) mass contours are shifted to larger values of mo and M1/2. This is in agreement 
with the previous observation that lower values of the seesaw scale lead to lighter sparticles 
for the same point in CMSSM parameter space (see discussion of Fig. [T|). As a consequence, 
the Higgs becomes lighter. To compensate for this downward shift in the SUSY spectrum 
one has to increase the parameters mo and/or M\ji- However, while a low type II scale 
now requires very large mo and/or M1/2, the resulting squark (and gluino) contours in the 
interesting range of niyfi are similar to those found in the CMSSM. This stems from the 
fact that m h o = 125 GeV requires again squark masses in the range of (at least) 5 TeV for 
Aq = OTeV and 2 TeV for A = — 5 TeV. This is not surprising since the Higgs mass is 



14 





FIG. 3. Squark mass, Higgs mass and Br(/i — > e^) contours in the plane (mo, -^1/2) for CMSSM 
plus seesaw type II for two values of the seesaw scale Mt' Mt = 10 9 GeV (top) and My = 10 14 GeV 
(bottom); as well as two values of Aq: Aq = OTeV (left) and Aq = — 5 TeV (right). The values 
of Br(/x — > e^f) are orders of magnitude below the expected experimental sensitivity in case of 
Mt = 10 9 GeV and, thus, are not shown. Note the change in scale compared to Fig. [2] (for further 
discussion see text). 



sensitive only to physical masses and mixings. However, as we will discuss below, there are 
some potentially interesting differences in the spectra due to the different RG running in 
the CMSSM and the SUSY type II seesaw. 

For the lowest value of Mt, where the spectrum distortions are larger, Br(/i — > e^y) is 
again negligible. Thus, an upper limit on Br(/i — > ej) provides an upper limit on Mt for 
any given value of the Higgs mass. A measurement of Br(/i — > cy) fixes a combination of 
A2 and Mt for fixed myfi. On the other hand, a lower limit on m h o provides a lower limit 
on a combination of m , Mi/ 2 and A for any fixed choice of Mt- Note that, contrarily to 
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FIG. 4. Dependence of several low-energy SUSY parameters on the couplings Ai j2 for the type 
II seesaw scenario. The results are shown for a specific point in the SUSY parameter space with 
m = 2.5 TeV, M 1/2 = 4.5 TeV, A Q = -5 TeV, tan/3 = 20 and fi > 0. Top left: Higgs mass as a 
function of A 2 for M T = 10 9 GeV (black), M T = 10 12 GeV (blue) and M T = 10 14 GeV (red). Top 
right: A 2 dependence of the stop masses (and their geometric average Ms), the Higssino mass 

parameter fi and the top-trilinear term At. Bottom left: stop mixing parameter Xt as function of 
A 2 - Bottom right: contours of the Higgs (black solid) and lightest stop (red dashed) masses in the 
(Ai, A 2 ) plane for M T = 10 9 GeV. The yellow regions corresponds to the CMS Higgs mass interval 
rriyfi = 125.3 ± 0.6 GeV. There are two CMS allowed contours in the lower right plot, since m h o 
first increases then decreases with A 2 for Mt = 10 9 GeV, compare to the figure in the upper left. 
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what happens in type I and III, in type II seesaw (with a single 15-plet pair) low-energy 
neutrino parameters essentially determine Yy in a way that large cancellations in the LFV 
soft masses are not possible. Moreover, when LFV in the soft masses is generated by Y^ 
only, the large value of sin #13 provided by the latest global analysis of neutrino oscillation 
data together with the present MEG bound on \x — > set an upper limit on the radiativ e 
r decays r — > /i(e)7 which is out of the reach of future experiments 



In the SUSY type II seesaw, the heavy triplet states T and T couple to the MSSM Higgs 
sector through the superpotential couplings A1.2 [see Eq. (TTUj) ]. We therefore expect these 
parameters to affect the Higgs mass to some extent. Obviously, since T and T are very 
heavy, the effect of Ai^ on the low-energy SUSY masses is indirect and originates from RG 
corrections induced on the SUSY parameters between Mqut and Mt- Consequently, these 
corrections are typically larger for smaller Mt- In Fig. H]we show the dependence of several 
parameters relevant for the computation of m h o as a function of Ai and A2, for a specific 
point of the SUSY parameter space (see caption). In the top-left panel we show a plot 
of niyfi versus A2 (and varying Ai from 0.1 to the maximum allowed by perturbativity) 
for Mt = 10 9 , 10 12 , 10 14 GeV. As expected, the impact of A2 on the Higgs mass is 
only significant for the case with Mt = 10 9 GeV. In the remaining two examples, a mild 
dependence on A2 is observed when the value of this parameter is very close to the Landau 
pole. From this plot one can also conclude that the effect of Ai on m h o is small, since the 
thickness of the lines (which reflects the variation of myp on Ai) is not too pronounced. In 
view of this, we will only comment on the A2-dependence of m^o for Mt = 10 9 GeV. 

The top-right panel of Fig. H] shows the variation of some relevant parameters with A2 
(see caption for more details). We first note that while \x and m ( - 12 (and consequently Ms) 
increase with increasing A2, \A t \ decreases (here A t is always negative since A < 0). This 
behaviour can be qualitatively understood by looking at the type II seesaw RGEs for the 
soft masses and trilinear terms. In particular, we notice that the RGE for A t contains a 
term proportional to | A2 1 2 ^4t which, at leading-log approximation, induces a correction to 
the top trilinear given by 

AA = -^^ Aoln {-MT-) ' (29) 

which is positive for A < 0. This explains why — A t decreases with A2 and, consequently, 
why X t decreases 0. The behaviour of \x can be traced taking into account that the Higgs 
soft masses rrii, receive a contribution which amounts to: 

11 d,u 

8vr 2 1 ' 1 V M T J 

Notice that, if not too small, the parameter Ai can act on rn 2 Hd as the top Yukawa coupling 
does on m 2 H bringing it to negative values at low-energies. In fact, we observe that for Ai 



Notice that in this case Xt ~ —At since // ~ — At and cot j3 = 0.05. 
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large, m 2 Hd is also negative at the EW scale. We recall from the EWSB symmetry breaking 
condition: 

" = Wl (31 > 

As A2 increases, m 2 Hu becomes more negative and m 2 Hd decreases, going from positive to 
negative values. This leads to an increasing of the value of /x with A2. 

Although not affected directly by A2, the stop masses mi (and, thus, the dynamical scale 
Ms) increases with that parameter mainly due to positive RG corrections in (nig 2c )33- The 
results for \X t /Ms\ as function of A2 are shown in the bottom- left panel of Fig. HJ Together 
with Eq. (127)1 . these results allow us to understand the behaviour of m h o with A2 shown in 
the top-left panel. In particular, we stress that for A2 — 0.5 we have \X t \ = a/6Ms, which 
corresponds to the "maximal mixing" scenario with maximised m h o (see the top- left panel). 
Finally, in the bottom-right panel, the contours of m^o and are shown in the Ai,2 plane 
for Mt = 10 9 GeV. The results confirm that while both m^o and depend reasonably 
strong on A2, their dependence on Ai is almost negligible. In particular, the effect of A2 on 
m h o can be much larger than its theoretical uncertainty. 

In Fig. [5] we show the results in the (mo, M1/2) plane for type III seesaw with Ms = 
10 13 , 10 14 GeV. For lower values of the 24-plet mass no solutions consistent with perturba- 
tivity exist. Since in type III all SUSY masses run strongly towards smaller values when 
Ms is lowered, already for Ms = 10 13 GeV the spectrum distortions with respect to the 
type I case are as large (or larger) as those found for type II with Mt = 10 9 GeV (compare 
Figs. |3]and [5]). As in type I and II, multi-TeV squarks (and gluinos) are required to explain 
m/jo ~ 125 GeV. Still, depending on Ms, the relations among sparticle masses are changed. 
It i s int eresting; to note that Br(/i — > ery) provides a particularly strong constraint for type 



III 107j. In case of Ms = 10 14 GeV (bottom plots in Fig. [5]), an improvement of Br(/x — > ej) 
to the level of ~ 10~ 13 (within the reach of MEG) would exclude the type III seesaw with 
degenerate 24-plets, R = 1 and a Higgs mass lying in the ATLAS and CMS range. In the 
particular case of Aq = —5 TeV (bottom-right plot) most of the CMS preferred region (in 
yellow) is already excluded by the constraint from MEG, which also excludes squark masses 
below ~ 7 TeV. 

We now turn to a discussion on differences found in the physical masses for the different 
seesaw setups. Since for type I the spectra are practically the same as in the CMSSM (which 
has been discussed at length in the literature) we focus on type II and III seesaws in the 
following. In Fig. [6] we show the allowed ranges for the lightest stop mass mi and the 
average squark mass, defined in Eq. (128]) . for seesaws of type II and III with a Higgs mass 
in the CMS range 125.3 ± 0.6 GeV. The red regions are for M^s = 10 14 GeV and the green 
ones for M T = 10 9 GeV and M E = 10 13 GeV. The left (right) plots are for A = OTeV 
(Aq = —5 TeV). Due to the CMSSM assumptions, stop and squark masses are tightly 
correlated, once the Higgs mass is fixed. It is interesting to note that once Aq is also set, 
the requirement that the Higgs mass falls into the CMS window leads to mass combinations 
which show a clear dependence on the seesaw scale. Especially noteworthy is the fact that 
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FIG. 5. Examples of squark mass, Higgs mass and Br(^i — > cy) contours in the plane (mo, M^) 
for CMSSM plus seesaw type III for two values of the seesaw scale M$s'- Ms = 10 13 GeV (top) 
and Mss = 10 14 GeV (bottom); as well as two values of Aq: Aq = OTeV (left) and Aq = — 5 TeV 
(right). For a discussion see text. 



no overlap between the regions with Ms = 10 14 GeV and Ms = 10 13 GeV exists in case of 
A = —5 TeV. Similar allowed mass ranges are obtained for seesaw type II. However, in 
this case we observe some overlap between the combinations of masses, even for the extreme 
cases of seesaw scales shown. It is nevertheless interesting that type III with a scale as low 
as 10 13 GeV does not allow squark and stop masses as large as type II does. For large values 
of Aq and fixed m h o, part of the parameter space is testable at the LHC with -^/s = 14 TeV. 
However, the allowed combinations of squark and stop mass for A = TeV are completely 
out of range of LHC14. 

Since in all our different setups large squark (and gluino) masses are required in order 
to explain a 125 GeV Higgs, the expectations are that no direct signals for SUSY will 
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FIG. 6. Allowed ranges of the lightest stop versus squark mass compatible with a 125.3 ± 0.6 GeV 
Higgs (CMS range) for type II (top) and type III (bottom) seesaws with Aq = TeV (left) and 
A = -5 TeV (right). The red (green) regions are for M T>S = 10 14 (M T = 10 9 GeV and M s = 10 13 ) 
GeV. 

be found in the near future. The LHC reach for i/i = 14 TeV and 300 (3000) fb" 1 has 
been recalculated very recently in 118| . The main conclusions of this study are that, via 



gluino/squark searches, LHC14 will be able to explore SUSY masses up to m- g ~ 3.2 TeV 
(3.6 TeV) for m~ q ~ m § and of m~ g ~ 1.8 TeV (2.3TeV) for m~ q > m g with 300 fb" 1 (3000 
fb" 1 ). Thus, for m^o ~ 125 GeV, only a small part of the allowed parameter space will 
be probed. However, future plans for the L HC envisage the possibility of ramping up the 



center-of-mass energy to a/s = 33 TeV |119| . With such a huge gain in energy, considerably 
larger regions of the parameter space allowed in our examples would become testable. 

To check this more quantitatively, we have calculated the cross sections for SUSY produc- 



tion at -y/i = 33 TeV for some representative points using the code Prospino 120 IB The 



5 The calculation of SUSY cross section at such large c.m.s. energy requires extrapolation of the measured 
PDFs and, therefore, is probably only a rough estimate. 
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m [TeV] M 1/2 [TeV] A [TeV] tan/3 sign(/x) M 5S [GeV] 



Point 1 3 3 20 + 10 14 

Point II 7 7 20 + 10 9 

TABLE I. Benchmark points with heavy squarks and gluino. Point I corresponds to a type I seesaw 
and point II to a type II seesaw. Both points have been chosen to give a Higgs mass of approximately 
m h o = 125 GeV. 



Particle Point I Point II 



X? 




1.35 


0.54 


d L , s L u L , 


CL 
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dR, SR UR, 


CR 
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bi 
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h 




5.9 


7.1 


h 




4.7 


5.0 


h 




5.6 


6.1 


9 




6.2 


2.7 



TABLE II. Some SUSY masses for the benchmark points given in Table |TJ All masses are given in 
TeV. For point I (point II) we find m h o ~ 125.6 (125.1) GeV. 



Production cross-section Point I Point II 





3.47 


2.08 


q q 


8.36 


0.60 


q q 


72.6 


9.59 


q g 


41.0 


793 


9 9 


3.49 


17000 



TABLE III. Most relevant production cross-sections for the benchmark points given in Tab le U All 
cross-sections are given in attobarns. These numbers have been computed with Prospino 



b lelfl . 



input parameters for two benchmark points lying inside the CMS Higgs mass range are given 
in Table HI We have chosen one point for type I seesaw (point I) and one for type II seesaw 
(point II), although for the SUSY production cross sections only squark and gluino masses 
are really important, of course. The corresponding SUSY spectra and some production cross 
sections are given in Tables ILT1 and ILTT| respectively. Point I has been deliberately chosen to 
give rriq ~ m g ~ 6 TeV, while point II leads to a heavier squark spectrum (around 7 TeV) 
but a lighter gluino. From Table III 1 1 one can see that point I (point II) would yield around 
~ 40 (~ 5300) squark/gluino events for an integrated luminosity of 300 fb~ x . These numbers 
are without any cuts and, therefore, should be taken as rough estimates. Nevertheless, they 
serve to illustrate how LHC33 would be able to cover most of the region of interest. This is 
also confirmed by Fig. [7] where we show the allowed regions in the (rriq, rrig) plane with m h o 
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m~(TeV) m~(TeV) 

FIG. 7. Allowed regions in the {mq,rrig) plane for seesaw type II (top) and III (bottom). In the 
left (right) panels Aq = OTeV (^4o = — 5TeV). We consider two extreme values of the seesaw scale 
in all cases as well as a fixed interval for the Higgs mass, namely m^o = 125.3 ± 0.6 GeV. For a 
discussion see text. 

in the CMS interval, and two extreme values of Mss, for type II seesaw (top) and type III 
seesaw (bottom). In the left (right) panel A = OTeV (A = — 5TeV). 

As before, we conclude that different seesaw models lead to distinct allowed combinations 
of masses. Still, in these plots, large overlaps between the regions for fixed A and different 
Mss are observed. Nevertheless, we find it especially encouraging that in type II and III 
seesaws gluino masses should be within the reach of LHC33 in almost all cases, for a Higgs 
mass in the CMS preferred window. The results also show that when Mt = 10 9 GeV 
(M T = 10 14 GeV), m g < 4.0 TeV (m s < 6.4 TeV) for A = OTeV. Instead, smaller values 
for the gluino mass are found if A < 0. The corresponding numbers for type III seesaw are 
m~ g < 5.1 TeV (m g < 5.4 TeV) for M 2 = 10 13 GeV (M s = 10 14 GeV) and A = TeV. These 
values should be compared with those of type I seesaw/pure-CMSSM where gluino masses 
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can be as large as m g < 7 TeV for Aq = TeV. 

A word of caution should be added to this discussion, owing to the fact that the upper 
limit on m g shown in Fig. [7] is very sensitive to the choice of the range for m^o. In particular, 
if the Higgs mass is as large as m^o = 128 GeV, which is currently not excluded, gluino masses 
up to 10 TeV and larger, would be allowed. Also, for small values of tan /3, say in the window 
tan /3 ~ (1 — 7), loop corrections to the Higgs mass are known to be small. This would again 
require much heavier stops and, therefore, much heavier gluinos to explain a m h a ~ 125 
GeV. 



V. CONCLUDING REMARKS 

In this work we have computed the mass of the lightest Higgs boson within the three tree- 
level realizations of SUSY seesaws and studied the main features of these models in light of 
the recent ATLAS and CMS results on Higgs mass searches. We have also complemented 
our analysis by considering the MEG bound on the LFV radiative decay u, — > ej. As in 
the pure CMSSM case, in SUSY seesaws a Higgs mass in the range (125 — 126) GeV (as 
preferred currently by CMS and ATLAS [3j, |4J) requires in all rather heavy SUSY 

spectrum. This is expected since m/jo is only sensitive to low-energy masses and mixings, 
and not to high-energy seesaw parameters (at least in a direct way). In other words, one 
can in principle find a different set of input parameters for each seesaw model leading to the 
same value of the Higgs mass. For this reason, a possible seesaw discrimination cannot rely 
on the Higgs mass data alone. Still, one expects to observe some differences in the physical 
low-energy SUSY spectrum. 

We have discussed squark, stop and gluino masses preferred by the current Higgs data in 
the different seesaw scenarios. While some small part of the parameter space allowed by a 
hefty Higgs will be tested at LHC14, most of our points are beyond the reach of th e next 

would 
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LHC run. However, a possible increase of the LHC energy to -y/i = 33 TeV 
make it possible to cover a large part of the parameter space allowed by the current Higgs 
data in our models. By considering some benchmark scenarios, we have also concluded 
that, in some cases, the allowed regions in the squark/stop and squark/gluino planes do 
not overlap when different values of the seesaw scale are considered or distinct seesaws are 
compared. Although this is not a general feature of the models under study, we believe this 
kind of analysis may be useful in the future to distinguish among seesaw setups and/or set 
limits on the input parameters of a particular model. Complementary information coming 
from the flavour sector, namely from rare decay searches, can also play a crucial role in 
the accomplishment of this task. In particular, upcoming data from MEG (and also from 
other LFV dedicated experiments) will certainly lead to further restrictions on the seesaw 
parameter space. 

We would like to mention that current data 0, 0] prefers an enhanced branching ratio 
for the di-photon final state; a obs /a SM = 1.54 ± 0.43 for CMS 1.9 ± 0.5 for ATLAS. With 
our heavy SUSY spectrum such an enhancement can not be explained. However, currently 
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this "discrepancy" is only of the order of (1 — 2) a and thus not significant. 

Finally, we would like to remark that, although at low-energies the CMSSM may not 
seem very different from its seesaw variants, the reconstruction of the initial conditions do 
drastically change from one case to the other. In view of this, one should reflect about the 
meaningfulness of fitting the CMSSM input parameters in a context where neutrino masses 
cannot be explained, as it happens to be in the MSSM. Low-energy measurements do result 
on different preferred regions for the input parameters when distinct models are considered. 
Obviously, this is not relevant for phenomenological studies at low energies, but it is surely 
crucial for studies addressing the dynamics behind SUSY breaking. 
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